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Abstract 

The width AT of the deconfinement transition region is shown to influence strongly 
the flow structure in the (Landau-) hydrodynamical expansion of a quark-gluon plas- 
ma. For a sharp first order transition (AT = 0) the mixed phase is rather long- 
lived, with a lifetime that has a maximum when the initial energy density is at the 
phase boundary between mixed and pure quark-gluon matter. For increasing AT, 
however, the lifetime decreases rapidly. Hadronic matter, however, remains long-lived 
as a consequence of the rapid change in the degrees of freedom in the transition region 
and the corresponding "softening" of the equation of state. 
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1 Introduction 



The hydrodynamical approach has been applied extensively to describe the dynamics of 
relativistic heavy-ion collisions ||, [3| , because it provides the only direct link between the 
equation of state of hot and dense nuclear matter and observable collective flow phenomena. 
Hydrodynamics is defined by local energy-momentum and (net) charge conservation, 

d fl T'» = Q, d^ = 0. (1) 

T^ u is the energy-momentum tensor, the (net) four-current of the ith conserved charge, 
IV = (iVf , N£, N£) for m conserved charges. Under the assumption of local equilibrium 
(the so-called "ideal fluid" approximation) the energy-momentum tensor and the m (net) 
charge currents assume the particularly simple form 0] 

T^ u = (e + p) u^u u —pg^ u , N? = run* 1 , i = l,...,m, (2) 

where e, p, are energy density, pressure, and (net) density of the zth conserved charge 
in the local rest frame of the fluid, w M = 7(1, v) is the four-velocity of the fluid (v is the 
three-velocity, 7 = (1 — v 2 ) -1 / 2 , = 1), and g^ v = diag(+, — , — , — ) is the metric tensor[j 
The system of m + 4 equations ([!]) is closed via eliminating one of the m + 5 unknowns in eqs. 
(0) by choosing an equation of state e.g. in the form p = p(e, n 1; ...,n m ). In the ideal fluid 
approximation, the (equilibrium) equation of state is the only input to the hydrodynamical 
equations of motion (|I|) that relates to properties of the matter under consideration and is 
thus able to influence the dynamical evolution of the system. The final results are uniquely 
determined once a particular initial condition and a decoupling ( "freeze-out" ) hypersurface 
are specified. 

In ultrarelativistic heavy-ion collisions (-y/s ~ 200 AGeV), the central region (in the 
center-of-momentum frame) is essentially free of (net) charges, due to the insufficient amount 
of baryon stopping H, and the charge conservation equations in (fj) do not have to be solved 
explicitly to describe the fluid evolution in this particular region of space-time. Moreover, 
the equation of state depends only on a single independent thermodynamical variable, for 
instance the temperature T. From the pressure as a function of T, p{T), one derives the 
entropy density s = dp/dT, and the energy density e = Ts — p. 

For this particular case, lattice QCD calculations provide fundamental information 
on the equation of state of strongly interacting matter. These calculations show a rapid 
increase of the entropy density as a function of temperature, s(T), in a narrow temperature 
interval AT ~ 10 MeV around a temperature T c ~ 150 MeV, indicating the presence of a 
transition from hadronic degrees of freedom at low temperatures to quark and gluon degrees 
of freedom at high temperatures. For a first order phase transition, the increase in s is 
actually a discontinuity, and AT is zero, while for a second or higher order transition the 
change of the entropy density is continuous, and AT is finite. Thusfar, however, lattice 



calculations have not converged on the order of the QCD phase transition H || [10], and we 
only know that AT < 0.1 T c . 



1 Our units are U 
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In this paper we discuss the sensitivity of the structure of hydro dynamic flow to this 
uncertainty in the equation of state. To this end we treat the width AT of the transition 
region as a free parameter and assume the following functional form for the entropy density 
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s(T) = c H T^- Unh " r 2 - T ^ + c Q ^ 1 + tM ' h '< r 2 - r °> /Ar ' . (3) 

Other functional forms are, of course, also possible and eq. (0) is chosen only for the sake 
of simplicity For T <C T c , the entropy density becomes that of a hadron gas, s(T) ~ chT 3 , 
while for T ^> T c , it assumes the value for a quark-gluon plasma, s(T) ~ cqT 3 . We denote 
the ratio of degrees of freedom in the quark-gluon and hadron phases by r = cq/ch- If 
the hadron phase is taken to be an ideal, massless pion gas, and the quark-gluon plasma to 
consist of (non-interacting, massless) u and d quarks and gluons, then r = 37/3. In the limit 
AT — > the equation of state (|[) features a first order phase transition between the two 
phases at the phase transition temperature T c and reduces to the MIT bag model equation of 
state [|l^]. In this paper we vary AT between and 0.1 T c , consistent with the uncertainties 
in present lattice QCD data. 

The first study of the sensitivity of hydrodynamic flow on AT was done in |TTJ . It was 



found that for AT < AT*, where AT* ~ 0.07 T c , matter is thermodynamically anomalous 
in the phase transition region and, consequently, that rarefaction shock waves occur in the 
one-dimensional expansion Jl3, [14| . On the other hand, for AT > AT*, matter is thermody- 
namically normal everywhere and simple rarefaction waves form the hydrodynamically stable 
solution. The main result of Ref. [|TTJ , however, was that, for the particular initial conditions 



chosen in that work, those structural differences had little effect on the final flow pattern. 
This led the authors of [|1T|] to the conclusion that observables, such as the distribution of 
outgoing particles, would hardly be affected by the uncertainty in AT. 

Renewed interest in the sensitivity of hydrodynamic flow on the equation of state has 



arisen as a consequence of the work of Refs. |L4], 15 , [L6| . It was demonstrated that a phase 



transition leads to a so-called ]T3| "softest point" in the equation of state, defined by a 
minimum in the ratio p/e. For the equation of state (0) with AT = 0, this point occurs at 
the phase boundary between mixed and pure quark-gluon matter (see also Section 2). It 
was argued |]15| that at this point the comparatively small pressure prevents a fast expansion 
and cooling of the system. This leads to a maximum in the lifetime of mixed phase matter as 
a function of the initial energy density. On this basis it was suggested |15j that the resulting 
enhanced emission of electromagnetic probes from this long-lived mixed phase could serve 
as a signature for the creation of the quark-gluon plasma in ultrarelativistic heavy-ion 
collisions. 



In Ref. fTql , the suppressed tendency for expansion at the softest point was shown to 
cause a minimum in the excitation function of the directed transverse collective flow in semi- 
peripheral collisions. It was proposed that this minimum would also be a unique signature 
for the phase transition. Another important observable sensitive to the lifetime of the plasma 
is pion interferometry [17]. In the standard "slowly burning log" scenario inspired by the 



small deflagration velocities associated with sharp transitions [18| one expects a dramatic 



difference between "out"- and "side" -correlations, which could serve as a signal of the phase 
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transition. 

In the present work, we investigate the hydrodynamic flow pattern for one-dimensional 



expansion in greater detail than in |ITJ . In particular, we focus on the question of the lifetime 
of different parts of the system as a function of the initial energy density eo and the width 
AT of the transition region. While one-dimensional hydrodynamics is not realistic for a 
detailed comparison with heavy-ion data, its simplicity allows us to study the influence of 
the equation of state on the hydrodynamic flow without additional effects originating from 
more complicated geometries. 

We confirm the observation made in [11 that, for initial energy densities well above or 
below the softest point of the equation of state, the temperature profiles in the thermody- 
namically normal case (AT > AT*) look indeed similar to those in the thermodynamically 
anomalous case (AT < AT*). However, we demonstrate that the temperature profiles are 
quite sensitive to a variation of AT for eo near the softest point of the equation of state. In 
particular, the lifetime of matter with temperatures near T c decreases rapidly for increasing 
AT. On the other hand, the lifetime of matter with temperatures well below T c remains 
large and may even grow for increasing AT. This is due to the "softening" of the equation 
of state by the change of degrees of freedom across the transition region. Thus, the de- 
tectability of the quark-gluon plasma via enhanced emission of electromagnetic probes may 
depend sensitively on AT, while other observables such as flow and correlation functions are 
expected to remain fairly insensitive to that uncertainty. 

This paper is organized as follows. In Section 2 we review properties of the equation 
of state (Q) in greater detail and establish the connection to thermodynamically normal or 
anomalous behaviour. Section 3 presents the expansion solutions as functions of eo and AT. 
In Section 4 we discuss the lifetimes of parts of the system having different temperatures. 
Section 5 concludes this work with a summary of our results. 



2 Equation of state and thermodynamical behaviour 



In this section we discuss the equation of state (||) and its thermodynamical properties. The 
entropy density at the phase transition temperature is s c = s(T c ) = cuij + l)^ 6 3 /2. With 
the definition p = (r + l)/(r — 1) and O(T) = (T — T C )/AT, one can write |l I ] 



-CO 



T 

T 



(l+p- 1 tanhe(T)) 



and subsequently obtains pressure and energy density from 

p(T) = [ T dT> s(T) , e(T) = T s(T) - p(T) 
Jo 

The square of the velocity of sound is 



dp _ s dT _ / 
de"Tds ~ [ 6 + 



T 



-i 



AT p + tanh0(T) cosh 2 Q(T) t 



(4) 



(5) 



(6) 
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It is obvious that for T C T c or T > T c (or, equivalently, |6(T)| ^> 1) the velocity of 
sound approaches the value for an ultrarelativistic ideal gas, c s = 1/ \/3. In the transition 
region, c s decreases and has a minimum. In the limit AT — > 0, c s approaches zero at 
T = T c and assumes the ideal gas value elsewhere. In Fig. 1 we show the energy density, 
the pressure, the entropy density, and the square of the velocity of sound as functions of 
temperature. Note that above T c , the energy density approaches the Stefan-Boltzmann 
value for a quark-gluon gas from above, while the pressure approaches it rather slowly from 
below. Both phenomena are observed in lattice data M, and are a consequence of general 
thermodynamical relationships |fl9f| . 

A way to classify the structure of one-dimensional hydrodynamical solutions [13, |l4]] is 
through the sign of 

d 2 p , 2c 2 (l-c 2 ) _ 



de 2 



+ 



e + p 



(An equivalent classifying quantity is the curvature of the Poisson adiabat [ p"T|j . 
equation of state @ we obtain 



For the 



3c 2 (l-3c 2 ) 



Matter is called thermodynamically normal if S is positive and thermodynamically anoma- 
lous if it is negative. It was shown in fll , |l3l [14| that in the first case the hydrodynamically 
stable solution for the one-dimensional expansion of semi-infinite matter is a simple rare- 
faction wave, while in the latter case it is a rarefaction shock wave. Depending on the value 
of AT, the equation of state (|3|) features thermodynamically normal as well as thermo- 
dynamically anomalous regions, as discussed below. The hydrodynamical solution for the 
one-dimensional expansion of semi-infinite matter may thus consist of a sequence of simple 
waves and shocks. For the one-dimensional expansion of finite matter, the expansion solu- 
tion is that for semi-infinite matter as long as the rarefaction waves have not reached the 
center of symmetry of the system (usually taken to be the center of the coordinate system, 
x = 0). Afterwards, the solution no longer consists only of simple waves and shocks but 
is more complicated (see also next section). An analytic solution is so far only known for 
thermodynamically normal matter with a constant velocity of sound [^0] . 

To see where matter described by ([3D becomes thermodynamically anomalous, it is ad- 
vantageous to plot thermodynamic quantities as a function of the energy density, since in 
the case AT = 0, the mapping T — > e(T) is not one-to-one. Moreover, the beam energy 
fixes the initial energy density, while the initial temperature is determined only indirectly 
through the equation of state. In Fig. 2 we show the entropy density, the ratio p/e, the 
velocity of sound (squared), and ETs as functions of e. Fig. 2 (b) serves to elucidate the 
notion of a "softest point" of an equation of state. For AT = 0, p/e = 1/3 = c 2 = const. 
in the hadronic phase. At the phase boundary between hadronic and mixed phase, e# = 
3T c s c /2(r + 1) = 0.1125 T c s c (for our choice of r = 37/3), the ratio p/e starts to decrease, 
since p = p c = T c s c /2(r + 1) = const, in the mixed phase while e increases. At the phase 
boundary between mixed and quark-gluon matter, eq = (4r — 1) T c s c /2(r + 1) = 1 .8125 T c s c 
(for r = 37/3), p/e assumes a minimum which is the so-called "softest point" 



15 . In the 
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quark-gluon phase, p/e increases again and reaches the asymptotic value 1/3 = c 2 for e — > oo. 
For increasing AT, the softest point is shifted towards smaller values of e. In addition, the 
equation of state becomes "less soft" at that point, i.e., the minimum value of p/e increases. 
On the other hand, there is a "softening" at energy densities around e#. 

Figs. 2 (c,d) show that, for AT = 0, the velocity of sound and S vanish in the mixed phase, 
tu < e < 6q. Although matter with a vanishing E is not thermodynamically anomalous 
in a rigorous sense, this nevertheless suffices to produce a rarefaction shock wave in the 
expansion |TJL |2T||. For AT being finite, but smaller than a critical value] 2 ] AT* ~ 0.07676 T c , 
there are always regions where £ becomes negative, i.e., where the equation of state is 
thermodynamically anomalous, and consequently, where rarefaction shock waves appear in 
the expansion. For AT larger than AT*, matter is thermodynamically normal everywhere 
and simple rarefaction waves form the hydrodynamically stable solution. Nevertheless, since 
the velocity of sound varies appreciably in the transition region, Fig. 2 (c), there will be 
modifications to a simple wave profile calculated with a constant velocity of sound. 

Note especially in Fig. 2 (d) that the region where matter becomes thermodynamically 
anomalous shrinks drastically even for tiny AT. To understand this we first expand the 
velocity of sound squared (|5|) for small AT/T C <C 1 in the vicinity of T = T c , 



Note that Q(T) = (T — T C )/AT is not necessarily small in this limit, i.e., one must not 
expand the hyperbolic tangents in this expression. In accord with Figs. 1 (d) and 2 (c), c 2 
is of the order AT/T C <C 1 near T c . Next, we determine the roots of eq. @ to lowest order 
in the small quantities AT/T C and T/T c — 1. This yields a quadratic equation for tanhO(T) 
with only one physical solution tanhO(Ts) = p — (p 2 + 3) 1 / 2 ~ —0.9173, corresponding 
to the upper boundary of the region where matter is thermodynamically anomalous. The 
lower boundary cannot be obtained in our approximation, because c 2 is already large in 
the respective region of energy densities, cf. Figs. 2 (c,d). Since tanh0(T s ) is of order 
unity, an expansion of the hyperbolic tangent in (|^) is indeed not possible. Note that 
T s /T c - 1 = (AT/T C ) Artanh[p - (p 2 + 3) 1 / 2 ] ~ —1.572 AT/T C . Thus, matter with tempe- 
rature T = T c is always thermodynamically normal for finite AT. Finally, we calculate the 
energy density corresponding to the upper boundary to lowest order, 



where e c = T c s c — p c = (2r + l)T c s c /2(r + 1) = 0.9625 T c s c . (This is the point where all 
curves in Fig. 1 (a) intersect.) This result is remarkable in several aspects. First, the fact 
that tanhB(Ts) is not small leads to a large shift of away from the point e c already to 
zeroth order in AT/T C . We checked numerically that first order corrections in AT/T C do not 
change the value for es appreciably before our approximations break down anyway. Second, 

2 This value was found numerically using eq. (^). 




AT/T C p + tanhO(T) 
T/T c 1 -tanh 2 9(T) ' 



(9) 



es - 




(10) 



6 



this shift is negative, i.e., away from eg which is the upper boundary in the case AT = 0. 
Thus, the size of the region where matter is thermodynamically anomalous does not change 
continuously in the limit AT — > 0, but increases abruptly when going from arbitrarily small 
AT to AT = 0. In this sense, the case AT = is very special in that it features an 
exceptionally large region where matter is thermodynamically anomalous. 



3 One— dimensional expansion solutions 

In this section the one-dimensional hydrodynamic expansion of a finite slab of matter with 
(homogeneously distributed) initial energy density eo is investigated with emphasis on the 
dependence of the flow structure on eo and AT. To solve the hydrodynamic equations (|l|) 
in (l+l)-dimensions we employ the relativistic Harten-Lax-van Leer-Einfeldt (HLLE) al- 
gorithm presented in |14|, |22| with a grid spacing of Ax = 0.025 R (2 R being the initial 
size of the slab) and a time step width of At = 0.99 Ax. This algorithm was shown |14 
to accurately reproduce hydrodynamic flow profiles for analytically solvable test problems 
where the equation of state has thermodynamically normal as well as anomalous regions. 
In particular, for the choice Ax = 0.025 R, the algorithm accurately reproduced the expan- 
sion solution for a finite system of size 2 R consisting of thermodynamically normal matter 
with a constant velocity of sound. (This solution was first provided by Landau in purely 
analytical form p0| , but can also be easily constructed semi-analytically via the method of 



characteristics p3[ .) 

In the following Figs. 3 - 7 we demonstrate how the flow structure changes systematically 
in the range e# < £o < 10 €q and for various AT. Each figure corresponds to a fixed 
eo, and shows the time evolution of temperature, center-of-momentum frame (CM) energy 
density, and the space-time structure of isotherms^ in three columns corresponding to AT = 
0, 0.01 T c , and 0.1 T c . 



3.1 eo = €h 

In Fig. 3 we show the expansion solution for an initial energy density eo = e# = 3T c s c /2(r + 
1) = 0.1125 T c s c . For AT = 0, Figs. 3 (a,d,g), this energy density corresponds to the bound- 
ary between hadronic and mixed phase matter. Thus, thermodynamically anomalous regions 
of the equation of state are not reached during the expansion, and the expansion proceeds 
at first as a simple rarefaction wave. After this simple wave overlaps with its symmetric 
counterpart coming from the (— x)-direction, the hydrodynamical solution is identical to 
that given by Landau [20|. 



For AT = 0.01 T c , Figs. 3 (b,e,h), the expansion seems to proceed rather similar to the 
previous case. Fig. 2 (d) shows, however, that around eo = tn the equation of state has 
a region where matter is thermodynamically anomalous. Thus, a rarefaction shock wave is 
bound to form in the expansion. Since the anomalous region is confined to a very narrow 

3 For AT = all states of mixed phase matter have the same temperature T c , cf. Fig. 1 (a). In this case, 
we define the (T = T c )-isotherm as the set of space-time points where T has just infinitesimally dropped 
below T c H. 



7 



range about e#, cf. Fig. 2 (d), the shock is weak and difficult to resolve numerically. The 
best evidence for its existence can be found in Fig. 3 (e), in that the slope of the rarefaction 
wave travelling into the slab is steeper than in the corresponding Figs. 3 (d,f) where no shock 
occurs. 

For AT = 0.1 T c , Figs. 3 (c,f,i), matter is thermodynamically normal everywhere, cf. Fig. 
2 (d). Consequently, the expansion initially proceeds as a simple rarefaction wave, as for 
AT = 0. For such a wave, the matter velocity with respect to the wave profile equals the 
velocity of sound 0, [H. Thus, the head of the simple wave travels with the velocity of 



sound c s into matter which is initially at rest. In the respective range of energy densities, 
c s is considerably smaller for AT = 0.1 T c than for AT = 0, cf. Fig. 2 (c). Therefore, as can 
be clearly seen in Figs. 3 (c,f), it takes much longer for the simple rarefaction wave to reach 
the center x = than in Figs. 3 (a,d). The delayed rarefaction can be also seen comparing 
Figs. 3 (g) and (i): it takes about 4R for the system at x = to cool below T = 0.7 T c and 
about 6 R to cool below T = 0.5 T c , while the corresponding values for AT = are 3 R and 
5R. 

The reason why there is no (T = T c )-isotherm in Fig. 3 (h) and no (T = T c )- and 
(T = 0.9 T c )-isotherm in Fig. 3 (i) can be inferred from Fig. 1 (a): as long as the energy 
density is below e c = (2r + l)T c s c /2(r + 1) = 0.9625 T c s c , the temperature corresponding 
to that energy density is smaller for larger AT. Thus, for AT = 0.01 T c , already the initial 
temperature is below T c , and for AT = 0.1 T c even below 0.9 T c , cf. Figs. 3 (b,c). Therefore, 
there simply are no corresponding isotherms in the space-time diagram. 

Note that the temperature and CM energy density profiles show a minimum at x = for 



late times t. As was explained in ||14|| , this is essentially a relativistic effect. At finite x, matter 
is moving, while due to symmetry it is at rest at x = 0. Moving matter, however, experiences 
relativistic time dilation and is thus less diluted and still hotter at a given time in the CM 
frame than matter at rest. Of course, the degree of dilation a fluid element experiences is 
determined by its velocity. This velocity is, in turn, determined by the "stiffness" of the 
equation of state, i.e., the ratio of pressure to energy density, since that ratio influences the 
acceleration of fluid elements. Thus, for the typical energy densities occurring in Fig. 3, 
matter velocities are smaller for larger values of AT, since the equation of state is softer, cf. 
Fig. 2 (b). This explains why the minimum at x = (or the time dilation effect, respectively) 
is less pronounced (and occurs at later times) in Figs. 3 (c,f) than for instance in Figs. 3 
(a,d). 



3.2 e H < e < e Q 

Fig. 4 shows the time evolution of temperature and CM energy density, and the isotherms 
in the space-time diagram for an initial energy density eo = 30T c s c /2(r + 1) = 1.125 T c s c , 
i.e., e H < e < e Q = (4r - l)T c s c /2(r + 1) = 1.8125 T c s c . For AT = 0, Figs. 4 (a,d,g), the 
initial state is thermodynamically anomalous mixed phase matter. Therefore, rarefaction 
proceeds via a shock wave, in front of which mixed phase matter is at rest. Hadronic matter 
emerges from the shock in an accelerated state and is subsequently diluted by a simple 
rarefaction wave fL4| . Note that the shock velocity is quite small, leading to a long lifetime 
of ~ 26 R for mixed phase matter at x = 0. We can understand the value for the lifetime 



S 



from the expression for the rarefaction shock velocity, eq. (70) of Ref. [H|, which reduces, 
for e H < e < e Q , to 

v sh * ~ - ^ . * "(2.6 eo/e^)- 1 • (11) 
3V3eo/eH 

Therefore, we expect for e /e# = 2(r + l)eo/3T c s c = 10 a lifetime tuf e = ~ 26 i?, in 

accord with Fig. 4 (g). 

Again, the expansion for AT = 0.01 T c looks similar, especially the temperature profiles, 
cf. Figs. 4 (a) and (b). The difference is that now the initial energy density is well above 
the region where matter is thermodynamically anomalous, cf. Fig. 2 (d). Therefore, the 
expansion has at first to proceed as a simple rarefaction wave, until the energy density has 
dropped sufficiently to reach that region. Then, a rarefaction shock wave provides further 
dilution. This can be most clearly seen in the CM energy density profiles, Fig. 4 (e). The 
reason why the temperature profiles, Fig. 4 (b), do not clearly exhibit this behaviour is that 
T is proportional to a fractional power of e, so that the range of high energy densities is 
not well resolved when plotted as a function of temperature. A clear indication for this 
behaviour, however, can be found comparing the space-time diagram Fig. 4 (g), where the 
(T = 0.9 T c )- and (T = T c )-isotherms coincide, since both temperatures occur at the position 
of the shock, with Fig. 4 (h), where they are well separated due to the simple rarefaction 
wave preceding the shock. Note that this has a drastic influence on the lifetime of matter 
with T = T c . Such matter is thermodynamically normal for finite AT, as discussed in the 
context of Fig. 2 (d), and consequently rapidly diluted by the simple rarefaction wave. The 
lifetime drops by more than a factor of 2 as compared to Fig. 4 (g). On the other hand, the 
lifetime of matter with T = 0.9 T c remains comparable to that in the case AT = 0, since 
this temperature occurs on the rarefaction shock front, cf. Fig. 4 (b), as was also the case in 
Fig. 4 (a). We will discuss this phenomenon in more detail in the next section. 

For AT = 0.1 T c , Figs. 4 (c,f,i), matter is thermodynamically normal everywhere and the 
hydrodynamically stable rarefaction solution is a simple wave (until it reaches the center 
x = and overlaps with its counterpart from the (— x)-region). Nevertheless, due to the 
variation of the velocity of sound in the phase transition region, cf. Fig. 2 (c), the shape of the 
simple wave exhibits a "bump" faintly resembling a smeared version of the rarefaction shock 
occurring in the previous two cases. This bump, however, is moving rapidly in positive re- 
direction, in contrast to the rarefaction shock in Fig. 4 (b), which is approximately stationary, 
or that in Fig. 4 (a), which is even slowly moving in (— x)-direction. This has the effect that 
the hot parts of the system (T > 0.9 T c ) cool faster than in the two previous cases. This 
can be understood with the help of Fig. 2 (b) which shows that for the respective energy 
densities the equation of state is "stiffer" for increasing AT, thus increasing the system's 
tendency to expand and cool. 

On the other hand, the cooler parts (T < 0.7 T c ) keep their temperature longer for 
AT = 0.1 T c , cf. also Figs. 4 (g,h,i). The reason is that the equation of state for small energy 
densities (low temperatures) is softer than in the other two cases, cf. Fig. 2 (b). Thus, the 
acceleration and consequently the matter velocities are smaller, and therefore the rarefaction 
and cooling take longer. We finally mention that, since now eo > e c , the initial temperature 
increases with AT, cf. Fig. 1 (a), in contrast to the previous case, Fig. 3. 
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3.3 The softest point eo ~ €q 

Fig. 5 shows the expansion solutions for e = 50T c s c /2(r + 1) = 1.875 T c s c . For AT = 0, 
Figs. 5 (a,d,g), this energy density is just above the softest point of the equation of state, i.e., 
the phase boundary between mixed and pure quark-gluon matter, €q = (4r — l)T c s c /2(r + 
1) = 1.8125 T c s c . Consequently, a simple rarefaction wave precedes the rarefaction shock 
wave travelling into hot quark-gluon matter flj]]. (The simple rarefaction wave is not visible 
in Fig. 5 (a) due to the large time step chosen.) This has the consequence that matter is 
pre-accelerated before it enters the shock wave. In conjunction with the equation of state 
becoming increasingly stiffer for smaller energy densities in the mixed phase, cf. Fig. 2 (b), 
this has the further consequence that the CM energy density develops a minimum at x = 
for given t, cf. Fig. 5 (d), just as in Fig. 3. We note that, in agreement with the results of 
|14| , the rarefaction shock is almost stationary for an eo this close to the softest point of the 



equation of state. This leads to a (local) maximum for the lifetime of matter with T = T c 
at x = 0. The particular value in Fig. 5 (g) can be again understood with the help of the 
approximate formula flTlD which gives^] tnf e ~ R/\v s h\ — 43 R. 

The temperature profiles for AT = 0.01 T c , Fig. 5 (b), look again similar to that of 
Fig. 5 (a), although the shock wave is no longer stationary, but at first advancing in in- 
direction. The reason is that the shock strength is smaller than in the previous case (due 
to the much narrower region where matter is thermodynamically anomalous, cf. Fig. 2 (d)), 
which is especially obvious from comparing Figs. 5 (d) and (e). Therefore, the velocity of the 
shock, which is originally directed inwards, is too small to compete with the pre-acceleration 
of matter in the simple rarefaction wave preceding the shock, resulting in a net outward 
motion of the shock. This mechanism is also responsible for the positive shock velocity 
observed for higher €q in the cases AT = and 0.01 T c , see Figs. 6, 7 below. For later times, 
however, just before the high energy density matter in the center is completely consumed, 
the shock reverses its direction, due to the fact that the energy density (or pressure) in the 
interior has now become too small to keep it driving outwards. (This effect is also visible in 
Fig. 5 (a).) 

The energy density profiles, Fig. 5 (e), do not show the minimum at x = typical for 
relativistic time dilation of moving matter. The reason is that in this case the equation of 
state is soft for the respective energy densities, cf. Fig. 2 (b). As explained above, this leads 
to small (i.e. non-relativistic) matter velocities and the absence of an appreciable amount of 
time dilation. 

The expansion for AT = 0.1 T c , Fig. 5 (c,f) resembles the corresponding one in Figs. 4 
(c,f). In particular, the "bump" structure caused by the variation of the velocity of sound is 
again clearly visible. As in Fig. 4, however, the expansion differs appreciably from the one 
for AT = and 0.01 T c in Figs. 5 (a,b,d,e). 

The relative time scales of cooling are similar to the previous case eo = 1.125 T c s c , cf. 
Figs. 4 (g,h,i) with Figs. 5 (g,h,i). Note again the sharp decrease in the lifetime of matter 
with T = T c for AT = 0.01 T c as compared to the case AT = 0, while matter with T = 0.9 T c 
remains long-lived. For AT = 0.1 T c , also the lifetime of matter with temperature T = 0.9 T c 



4 We may safely neglect the small pre-acceleration of matter due to the simple rarefaction wave preceding 
the shock in the quark-gluon phase. 
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is reduced, while matter that is even cooler still lives longer. 



3.4 e ~2e Q 



Fig. 6 shows the expansion for an initial energy density eo = 100 T c s c /2(r + 1) = 3.75 T c s 
2€q. This case is of interest since for initial energy densities of this order of magnitude a 
rather peculiar phenomenon occurs in the expansion for an equation of state with AT = 0. 
One observes in Figs. 6 (a,d) that the rarefaction shock wave is now driven outwards by the 
pre-acceleration from the simple rarefaction wave, as was already discussed above. Matter 
behind the shock experiences cooling typical for the Landau expansion with a sufficiently 
stiff equation of state, i.e., slower matter near the center x = cools more rapidly than faster 
matter at larger x = 0. This, however, has the consequence that as soon as the (comoving) 
energy density drops below e#, or equivalently, the temperature drops below T c , a second 
rarefaction shock wave develops that connects hadronic matter in the center with mixed 
phase matter behind the first shock. 

The occurrence of the second shock could have been expected, since matter changes 
from thermodynamically anomalous in the mixed phase to thermodynamically normal in the 
hadron phase. To our knowledge there is now, however, no general mathematical theorem 
rigorously proving the existence of this second shock, as was the case in the expansion of 
semi-infinite matter studied in fH|. The shock has an interesting further consequence: just 



as the first shock expells hadronic matter outwards, the second accelerates matter inwards. 
Therefore, there is an accumulation of matter in the center which is consequently reheated 
instead of cooled. This reheating of the center can be clearly seen in Fig. 6 (g), where the 
triangle- shaped isotherm at small x and t ~ 70 i? is a re-occurrence of the (T = 0.9 T c )- 
isotherm caused by reheating. In addition, also the (T = 0.7T c )-isotherm shows that the 
central region remains hot over a much longer period of time than for the expansion of matter 
without a phase transition. 

The reheating is considerably smaller for AT = 0.01 T c , cf. Figs. 6 (b,e,h), and basically 
only discernible from the shape of the (T = 0.7T c )-isotherm. The reason is, of course, that 
the region in the equation of state, where matter is thermodynamically anomalous, is much 
narrower, cf. Fig. 2 (d), so that the strength of the second rarefaction shock is smaller and 
consequently the acceleration of matter towards the center much weaker. For AT = 0.1 T c 
there is no such an effect at all, and the fluid evolution is rather different from the cases 
AT = and 0.01 T,. 



3.5 e > e Q 

Finally, we show in Fig. 7 the expansion for eo = 500T c s c /2(r + 1) = 18.75 T c s c ^> €q. 
This case is similar to the one studied in Fig. 13 in Ref. UTT] , although we followed the fluid 
evolution over a much longer time. Here, the acceleration from the simple wave travelling 
into the quark-gluon phase is rather large, leading to a shock wave in Figs. 7 (a,b,d,e) 
that is rapidly moving outwards. Note that now also the temperature profiles for AT = 
0.1 T c , Fig. 7 (c), resemble those in the case AT = 0, Fig. 7 (a) (with the exception of the 
immediate vicinity of the rarefaction shock). The main reason is that the fluid motion is now 
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dominated by the conversion of the large initial energy density into collective motion, so that 
the influence of the phase transition region on the structure of the hydrodynamic solution 
becomes less important. Moreover, as in the cases considered previously, the variation of the 
velocity of sound in the phase transition region causes a "bump" structure that resembles a 
smeared shock. 

The similarity of the temperature profiles in Figs. 7 (a) and (c) led the authors of JTIj to 
the conclusion that final state observables will not be much influenced by the precise value 
of AT. They reached the same conclusion studying the expansion with an initial entropy 
density so = 0.3 s c (Fig. 11 in ||11|| ) which corresponds, for AT = 0, to an initial energy 
density eo = 0.2625 T c s c . For this case, the expansion is similar to the one studied here 
in Fig. 3. That figure shows that the influence of AT on the expansion is negligible also 



for such small initial energy densities, confirming the conclusions of |Tl| . As we have seen 



however, the profiles are rather sensitive to the choice of AT for initial energy densities near 
the softest point of the equation of state. In particular, the reheating phenomenon of Figs. 6 
(a,d,g) does not occur at all for AT = 0.1 T c , cf. Figs. 6 (c,f,i). We discuss the effect on the 
lifetime of matter with different temperatures in the next section. 



4 Lifetime of parts of the system having different tem- 
peratures 

In this section we discuss the lifetimes of parts of the system having different temperatures. 
The "lifetime" is defined as the time it takes the system to cool below a certain temperature 
T* at x = 0. More precisely, the "lifetime" is the intercept of the isotherm t(x; T*) in the 
(t — x)-plane with the t-axis, i.e., tnf e = t(0;T*). We choose the point x = to measure 
this time, since the motion of fluid cells at finite x causes a time dilation which leads to a 
slower cooling when measured in the global CM frame of the system. In case there is more 
than one intercept, for instance when there is reheating of the center, Fig. 6 (g), we choose 
the first. 

We investigate the lifetime in two different geometries. In subsection 4.1 we study longi- 
tudinal expansion for fixed initial "radius" R of the system. In subsection 4.2 we investigate 
a more realistic geometry where the initial longitudinal size R takes into account the Lorentz 
contraction of the colliding nuclei. Thus, R varies with beam energy, and hence with the 
initial energy density eo- 



4.1 Fixed R 

In Fig. 8 we show the lifetime of parts of the system having temperature T = 0.5 T c , 

0.7T C , 0.9 T c , T c , 1.1 T c , and 1.3 T c for the three different AT = 0, 0.01 T c , 0.1 T c , and for 

an ideal gas equation of state without phase transition and constant velocity of sound 

c s = 1/a/3 as a function of the initial energy density eo- We measure the lifetime in units 

of R. For the ideal gas calculation, Fig. 8 (d), there is of course no phase transition, and 

the temperature scale set by T c is arbitrary. The results of Fig. 8 (d) were obtained using 
T /T c =(4e /3T cSc ) 1 /4. 
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c 



One observes that for initial energy densities smaller than ~ e c , the lifetime of matter 
with temperatures T < 0.9 T c is essentially independent of AT. The only exception is that 
for AT = and 0.01 T c the lifetime of matter having T = 0.7 T c and T = 0.9 T c is equal, 
while that is not the case for AT = 0.1 T c . The reason is obvious from Figs. 4, 5: in the first 
two cases the rarefaction shock wave, once it reaches the center, abruptly cools matter to 
about 0.5 T c , thus terminating the existence of matter having T = 0.9 T c as well as T = 0.7 T ( 
simultaneously. This can be clearly seen in the space-time diagrams Figs. 4, 5 (g) and (h). In 
the case AT = 0.1 T c , rarefaction proceeds via simple waves which are continuous solutions 
and gradually cool the system, and thus matter with T = 0.7 T c survives longer than hotter 
matter with T = 0.9 T c . 

The second observation is that the lifetime of "very hot" parts of the system, i.e., matter 
having temperature T = 1.1 T c and 1.3 T c , is also comparable for all three values of AT. 
These high temperatures occur, of course, only for sufficiently high initial energy densities. 
The reason for the similarity in the lifetimes is that a variation of AT has little effect on the 
equation of state for high energy densities (or temperatures), cf. Fig. 1. Consequently, also 
the fluid evolution must be similar for these high temperatures . 

For initial energy densities near the softest point of the equation of state and for temper- 
atures around T c , however, AT has an appreciable influence on the lifetime. For instance, 
one clearly observes in Fig. 8 (a) that the lifetime of matter with T = 0.9 T c and T = T C has a 



maximum near the softest point, cf. Fig. 2 (b). This effect was discussed in detail in Ref. [|14 



and in the context of Fig. 5 above. It has been explained by the existence of a minimum in 
the velocity of the rarefaction wave travelling into the system, leading to the "slow burning" 
of hot matter also observable in Figs. 4, 5 (a,d). An alternative explanation was given in 
[D|, namely that the small p/e-ratio at this point reduces the system's tendency to expand. 

However, Figs. 8 (b,c) show that this maximum in the lifetime is rapidly washed out 
and eventually vanishes for increasing AT, although the p/e-ratio still exhibits a minimum, 
cf. Fig. 2 (b). This vanishing of the maximum in the lifetime of matter with a particular 
temperature is essentially explained by a gradual reduction of the lifetime when increasing 
AT: first, the lifetime of matter with T = T c is reduced when increasing AT from to 
0.01 T c , cf. Figs. 8 (a) and (b). Then, that of matter with T = 0.9 T c is reduced when further 
increasing AT from 0.01 T c to 0.1 T c , cf. Figs. 8 (b) and (c). The lifetime of cooler matter 
with T < 0.7 T c is not appreciably influenced around e = e<g. 

In order to understand this reduction of the lifetime for increasing AT, we show in Fig. 
9 the time evolution of the temperature at x = for various values of AT and an initial 
energy density e = 50T c s c /2(r + 1) = 1.875 T c s c . This value is close to the softest point 
of the equation of state, i.e., where the above mentioned effects on the lifetime are most 
drastic. In this representation, the rather long lifetime of matter with temperature T = T c 
(tufe ~ R/\v s h\ — 43 R) for a sharp first order phase transition can be seen in a particularly 
clear way. This temporal temperature profile changes smoothly with increasing AT. The 
drastic effect on the lifetime of matter with T = T c and T = 0.9 T c observed in Fig. 8 is 
readily explained via Fig. 9 as follows. For AT = the temperature stays constant (and 
equal to T c ) over a rather long period of time ~ 43 R, due to the small velocity of the 
rarefaction shock which consumes mixed phase matter and consequently takes a long time 
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to reach the center x = 0, cf. Fig. 5 (a). After that, the temperature drops sharply, cf. also 
Fig. 5 (g). 

For small finite AT = 0.01 T c , however, the temperature is already slowly decreasing with 
time, cf. also Figs. 5 (b,e). This leads to a much shorter lifetime of matter with the particular 
temperature T = T c , although the system remains obviously rather hot, T ~ T c , for a time 
span ~ 37 i? that is only slightly smaller than in the previous case. This also accounts for 
the long lifetime of matter with T = 0.9 T c as observed in Fig. 8 (b). After that time, the 
cooling is again rather drastic due to the rarefaction shock which has then reached the center 
x = 0. For increasing AT, the cooling of hot matter becomes faster, which explains why the 
lifetime of matter with T = 0.9 T c is eventually also appreciably reduced, cf. Fig. 8 (c). 

One also notes that the lifetime of matter with T = 0.7 T c is not much influenced when 
increasing AT, in accord with our above observation in Fig. 8. On the other hand, mat- 
ter with even lower temperatures may survive longer. In general, the more rapid cooling 
of matter with higher temperatures as well as the delayed cooling of matter with smaller 
temperatures is well explained by the increasing "stiffness" of the equation of state at higher 
energy densities and its "softening" at smaller energy densities for increasing AT, cf. Fig. 2 
(b). 

An important conclusion to be drawn from Figs. 8 and 9 is that for increasing AT 
the maximum in the lifetime of matter with temperature around T c vanishes and that, for 
AT =0.1 T c , there is no longer any irregularity that can be traced to the rapid variation of 
thermodynamical quantities in the phase transition region. In fact, in that case the shapes 
of the lifetime curves resemble those for the expansion of an ideal gas, cf. Figs. 8 (c) and (d). 
However, the existence of a transition from hadronic to quark and gluon degrees of freedom, 
the softening of the equation of state induced by that transition, and the resulting delayed 
expansion causes the lifetimes to be larger by an overall factor ~ 10, cf. also Fig. 9. We 
will discuss possible implications for proposed quark-gluon-plasma signatures in the next 
section. 

The curves in Fig. 8 (d) could also be calculated analytically employing Landau's solution 
20| which gives 

dx(T,a) 



t 



x=0 



dT 



;i2) 



x=0 



where 

TR rHTo/T] 



x (T,a) = r ° du e (1+/3)u I (f3Ju* - c 2 s aA . (13) 



t\ x=0 = - -r Io (HTa/T)) - / du e 2u I (u) . (14) 



Here (3 = (1 — c 2 )/2c 2 = 1 for c 2 s = 1/3 and a = tanhf is the fluid rapidity. At x = we 
also have a = by symmetry and we obtain 

To understand the qualitative behaviour of the curves in Fig. 8 (d), it is sufficient to note 
that the integral is always dominated by contributions from the upper boundary. For To 3> T 
one may thus use the asymptotic approximation Iq(u) ~ e u in the integrand and obtains 

t\ x=0 = t life ^RT^Re 3 /A (15) 
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for fixed T. Thus, tuf e / R behaves (at least asymptotically) ~ ej^ 4 independent of T/T c . 
This behaviour is in accord with the one observed in Fig. 8 (d). Curves for different T/T c 
are simply shifted due to different constants of proportionality. Note that an increase of the 
lifetime proportional to e^ 4 is also expected in scaling hydrodynamics |L4 . 



4.2 Variable R(e ) 

In Fig. 8 R was assumed fixed for all initial energy densities. Now we allow it to vary 
with eo to take into account the Lorentz contraction of the colliding nuclei. To specify this 
relationship quantitatively, we consider two models for the initial stage of central A + A— 
collisions. 

The first is the scenario employed in ||15|| . It is assumed that half of the available CM 
energy, \fs = 2 AMnJcm, is deposited in a cylinder of height IRaIicm and radius Ra- 
Thus, e = AM N ^£ M /2nR\. The initial "radius" of the system is therefore 

p Ra Ra ( AM N \ 1/2 

1CM (e /T c s c y/i\4n(r + l)p c R A ) ' 1 ] 

With p c = tt 2 T c 4 /30, T c = 170 MeV 0, and R A = 1.12 A 1 / 3 fm, this leads to R/R A ^ 
0.334 (e /T c s c ) -1 / 2 . Fig. 10 shows the results of Fig. 8 with the tuf e / R-scaie converted 
according to this formula. One observes that the lifetime in units of the radius Ra of the 
nucleus does not grow as rapidly with the initial energy density as in Fig. 8, due to the fact 
that the system's initial size is shrinking with the CM-frame Lorentz-gamma (leading to the 
factor €q 1//2 in (|i~6D ). This suppression of the lifetime at high e has the interesting effect that 
existing maxima in the lifetime of different temperatures as observed in Fig. 8 are enhanced 
in Fig. 10. In the case AT = 0.1 T c , even new maxima are created for T = 0.7 and 0.9 T c . 
The reduction of the lifetimes of matter with temperatures near T c for increasing AT is, of 
course, not influenced by the rescaling of the tuj e / R-axis. 

On the other hand, in the ideal g clS CclSC, Fig. 10 (d), the rescaling of the tuf e / R-ax.is 
creates a minimum. This can be understood as follows. For large eo we may employ the 
asymptotic formula ( |i~5"D and note that R ~ e l ' 2 according to (|i~6D, which results in tuf e ~ 
e^ 4 , i.e., a weak increase with eo as observed in Fig. 10 (d). For small eo, we employ ( |T4] ) 
with T ~ T , which results in tnf e ~ R ~ e , i.e., a decrease with eo- The two different 
forms give rise to the observed minimum at a certain value of eo which increases for increasing 
T/T c , since the asymptotic behaviour T ^> T is then realized only for correspondingly larger 
values of eo ~ Tq . The qualitative difference in the behaviour seen in Fig. 10 (c) as compared 
to Fig. 10 (d) contrasts sharply the similarity of Figs. 8 (c) and (d). The maxima in Fig. 
10 (c) and the minima in Fig. 10 (d) are, however, not prominent enough to be of relevance 
for quark-gluon-plasma signatures. On the other hand, the huge quantitative difference in 
the lifetimes (a factor of ~ 10), which is, of course, still present, could have consequences 
for observables (see next section). 
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This value is slightly larger than the one cited in the Introduction, but still consistent with lattice QCD 



calculations, cf. the discussion in M. We choose this value mainly to maintain consistency with [fL6, 24 
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For the second scenario, we employ the one-dimensional shock model UM. Since we 



presently consider baryon-free matter only, we model the "ground state" of the colliding 
nuclei as a (massless) pion gas at T ~ 143 MeV, i.e., with pressure pi = p c /2. The energy 
density e of the hot, compressed state is determined by solving the integral conservation 
laws for energy and momentum across the shock Q. This gives also values for the shock 
velocity v' sh in the rest frame of incoming matter 



v sh = — , (17) 

V e - Cj ti+po 



1/2 



and the CM velocity vqm of matter required to reach that particular energy density 

1/2 

eo + Pi £i + Po 



Po-Pi e - 6j . . , 

vcm = ; ; • (18) 



The CM time when the complete "nucleus" is compressed (or the shock front has traversed 
the incoming nucleus, respectively) is calculated as [[13], Q 

(19) 



t£_ = (1 - v' sh VCM)lCM 

R A v' sh 



This time is (up to a factor c = 1) identical to the initial "radius" R of the system. With 
equations ( |P7D and ( |i~8"|) one can thus relate R to e , similar as in (|T6|). One has, however, 
to consider that a single shock is not the hydrodynamically stable solution in thermody- 



namically anomalous matter pi |. The correct solution is a sequence of single shocks and 
simple compression waves and, as a consequence, the initial energy density of the system is 
no longer constant over its extension 2R. For our purpose of obtaining a simple estimate 
of R as a function of eo we will neglect these details and use a linear interpolation in the 
respective region of eo-values. The resulting initial size R of the system as a function of eo is 
shown in Fig. 11. With this R(eo) we then convert the tuf e / R-scale in Fig. 8. The results are 
depicted in Fig. 12 and are qualitatively similar to Fig. 10. The rapid decrease of R above 
eo — eg for AT = and 0.01 T c strongly enhances the maximum in the lifetime of matter 
with temperatures near T c and even creates a maximum in the lifetime of matter with cooler 
temperature, cf. Figs. 12 (a,b). The maxima observed in Fig. 10 (c) are now also more pro- 
nounced, cf. Fig. 12 (c), while Fig. 12 (d) still exhibits minima as in the corresponding Fig. 
10 (d). The one-dimensional shock model thus proves to be the most optimistic scenario 
for quark-gluon-plasma signatures that require the existence of a distinguished maximum 
in the lifetime of hot matter at the softest point of the equation of state, like the proposed 



enhanced emission of electromagnetic probes ITS 



We conclude that the lifetime of matter with temperatures near T c is rather sensitive to 
the choice of AT and decreases rapidly for increasing AT for initial energy densities around 
the softest point of the equation of state. On the other hand, regions of lower temperature 
(T < 0.7 T c , cf. Fig. 9) are less sensitive to variations in AT and even survive longer with 
increasing AT. Both phenomena have been explained in terms of the change in the "stiffness" 
of the equation of state when varying AT. For a fixed initial system size (independent of eo), 
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increasing AT eliminates the maxima in the lifetime near the softest point of the equation of 
state (see Fig. 8 (c)). On the other hand, if one accounts for the Lorentz contraction of the 
initial size of the system, then such maxima in the lifetime persist even for finite AT. This is 
in contrast to the case without a phase transition, where the lifetimes rather exhibit minima 
for the corresponding values of eo (Figs. 10, 12 (d)). The lifetimes are also much longer 
(a factor of ~ 10) in case a rapid change in the degrees of freedom occurs. Experimental 
implications will be discussed in the next section. 



5 Summary and Conclusions 

In this work we have studied the one-dimensional hydrodynamic expansion of a finite system 
with an equation of state featuring a transition between (baryon-free) hadronic matter (with 
3 massless degrees of freedom) and a quark-gluon plasma (with 37 massless degrees of 
freedom). The system's evolution was studied as a function of its initial energy density eo 
(which was assumed to be homogeneous) and the width AT of the phase transition region 
in the equation of state. For AT < AT*, where AT* ~ 0.07676 T c , the equation of state 
has regions where matter becomes thermodynamically anomalous. Thus, rarefaction shock 
waves appear in the one-dimensional expansion, in contrast to the case AT > AT*, where 
matter is thermodynamically normal for all energy densities and (overlapping) simple waves 
form the hydrodynamically stable expansion solution. 

Temperature profiles are not very sensitive to AT for initial energy densities well above 
and below the softest point of the equation of state, in agreement with the results of [ |TT 



Variations of AT have, however, an appreciable influence on the temperature profiles for 
eo near the softest point. In particular, rarefaction shocks are much weaker for finite AT 
than in the case AT = and transform into "bumps" resembling smeared shocks in the 
thermodynamically normal case AT > AT*. The shock strength is much smaller because, 
as shown via eq. ([T0|), the region where matter is thermodynamically anomalous shrinks 
abruptly when AT becomes finite. The "bump" structures in the thermodynamically normal 
case appear due to the variation of the velocity of sound in the transition region. In general, 
the equation of state at high energy densities becomes increasingly stiffer for increasing AT, 
which leads to more rapid expansion and cooling of the hot parts of the system. On the 
other hand, it becomes softer at smaller energy densities and thus the cooler parts of the 
system expand less rapidly and keep their temperature longer. 

The above results were obtained assuming one-dimensional geometry and neglecting the 
finite (net) baryon density in the stopping regime. Nevertheless, we consider possible exper- 
imental consequences that may apply also in more realistic situations. The most important 
consequence from the experimental point of view is that the lifetimes of matter with temper- 
atures near T c decrease rapidly with increasing AT, if eo is in the vicinity of the softest point 
of the equation of state. This will have an effect on the electromagnetic signature proposed 
in [|15|]. For instance, the rate (per unit volume) of direct thermal photons with energy E 
emitted from quark-gluon matter (at rest) is given by 
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For a fixed photon energy of E = 1 GeV, a reduction of the temperature by ~ 10 — 20%, 
as seen in Fig. 9 when AT increases from to 0.1 T c , would reduce the number of photons 
by a factor of ~ 2 — 4 (for T c = 170 MeV, as = 0.4). In addition, the rate integrated over 
space-time will be smaller by an even larger factor, because the space-time volume of hot 
matter also shrinks for increasing AT, cf. Figs. 3-7 (g,h,i). This makes a direct thermal 
electromagnetic signal of the quark-gluon plasma much more difficult to disentangle from 
background sources. On the other hand, after taking into account the expected Lorentz 
contraction of the initial size of the system, all lifetimes assume a maximum as a function 
of eo at the softest point of the equation of state, cf. Figs. 11, 12, independent of AT. The 
detection of this feature is, however, complicated by the broadening of this maximum with 
increasing AT. 

The generic feature that remains independent of AT is that the lifetimes of matter of 
any temperature are larger by a factor of ~ 10 as compared to the ideal g In this 

respect, the interferometric signal of such long lifetimes associated with a phase transition 
to the quark-gluon plasma W7\ appears to be less sensitive to AT than hard electromagnetic 



probes. 

Finally we comment on the signature proposed in Ref . JIB] , where it was shown that the 
excitation function of the directed transverse flow has a minimum generated by the softest 
point of the equation of state, or in other words, the minimum in the function p/e. Such 
a minimum still exists even for finite AT, cf. Fig. 2 (b). However, the intrinsic dynamical 
mechanism for the decrease of the transverse directed flow was the inertness of the mixed 
phase that prevented it to expand and deflect spectator matter. This inertness is weakened 
at finite AT and hot matter in the phase transition region expands and cools much more 
rapidly than in the case of a strong first order phase transition, although still much slower 
than if there were no transition at all, cf. Figs. 8 (c) and (d). Further and more detailed, 
three-dimensional studies |2H] are required before a definite conclusion of the effect on the 



transverse directed flow can be drawn. 
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Figure Captions: 



Fig. 1: (a) the energy density (in units of T c s c ) divided by (T/T c ) 4 , (b) three times the 
pressure (in units of T c s c ) divided by (T/T c ) 4 , (c) the entropy density (in units of s c ) divided 
by (T/T c ) 3 , and (d) the square of the velocity of sound as functions of T/T c . The solid line 
with the discontinuity at T = T c in (a,c) corresponds to AT = 0. In (b), this curve has a kink 
at T = T c , and in (d) it is the flat line at c 2 s — 1/3 (the point T = T c is an exception, there 
c 2 s = 0). The other curves correspond to AT = 0.01 T c (dotted), AT = 0.05 T c (dashed), 
AT = 0.075 T c (dash-dotted), and AT = 0.1 T c (solid). 

Fig. 2: (a) entropy density (in units of s c ), divided by (T/T c ) 3 , (b) p/e, (c) the veloc- 
ity of sound squared, and (d) STs as functions of energy density (in units of T c s c ). The 
different curves correspond to those in Fig. I. 

Fig. 3: (a,b,c) temperature profiles and (d,e,f) center-of-momentum frame energy den- 
sity profiles for times t = 0.75 nXR, n = 0,1,..., 7 and an initial energy density of e = 
3T c s c /2(r + 1) = 0.1125 T c s c = e H . The profiles are alternatingly shown as full and dotted 
lines in order to better distinguish them. (g,h,i) show isotherms in the space-time diagram. 
The isotherms are labelled with the corresponding temperatures in units of T c . Figs. (a,d,g) 
are for AT = 0, (b,e,h) for AT = 0.01 T c , and (c,f,i) for AT = 0.1 T c . 

Fig. 4: Same as in Fig. 3, for e = 1.125 T c s c , i.e., e H < e < e Q = (4r — l)T c s c /2(r + 1). 
Profiles in (a-f) are for times t = AnXR, n — 0, 1, 7. 

Fig. 5: Same as in Fig. 3, for e = 1.875 T c s c ~ eq. Profiles in (a-f) are for times 
t = 6n\R, n = 0, 1,...,8. 

Fig. 6: Same as in Fig. 3, for e = 3.75 T c s c ~ 2cq. Profiles in (a-f) are for times 
t = WnXR, n = 0, 1,...,7. 

Fig. 7: Same as in Fig. 3, for e = 18.75 T c s c ^> eq. Profiles in (a-f) are for times 
t = 12.5nXR, n = 0, 1,...,8. 

Fig. 8: The lifetime of differently hot parts of the system as a function of the initial energy 
density for (a) AT = 0, (b) AT = 0.01 T c , (c) AT = 0.1 T c , and (d) an ideal gas equation of 
state with c 2 s = 1/3. Full line and circles: lifetime of matter having temperature T = 0.5 T c , 
dotted line and squares: T = 0.7 T c , dashed line and diamonds: T = 0.9 T c , long dashed line 
and triangles pointing upwards: T = T c , dash-dotted line and triangles pointing downwards: 
T = 1.1 T c , full line and stars: T = 1.3 T c . Note the change in the scale of tuf e in (d). 

Fig. 9: Temperature (in units of T c ) at x = as a function of time (in units of R) for 
e = 1.875 T c s c ~ £q. Lines correspond to AT-values as in Figs. 1, 2. The dashed line 
corresponds to a calculation with an ideal gas equation of state, cj. = 1/3. 



21 



Fig. 10: The same as in Fig. 8, except that tuf e /R is converted to tnf e /RA according 
to eq. (|l6f) , where Ra is the radius of a nucleus at rest. 



Fig. 11: The system's longitudinal "radius" R in units of Ra as a function of the ini- 
tial energy density in the one-dimensional shock model. Full line: AT = 0, dotted line: 
AT = 0.01 T c , dashed line: AT = 0.1 T c , dash-dotted line: ideal gas equation of state. The 
straight parts of the full and dotted lines are a linear interpolation between the points where 
the hydrodynamically stable compression solution ceases to be a single shock due to the 
thermodynamically anomalous behaviour of the equation of state in the phase transition 
region. 

Fig. 12: The same as in Fig. 8, except that tuf e /R is converted into tnf e / Ra according 
to Fig. 11. 
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